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ABSTRACT 

The Gamma- Ray Burst (GRB) 021211 detected by the High Energy Transient Ex- 
plorer (HETE) II had a simple light-curve in the x-ray and gamma-ray energy bands 
containing one peak and little temporal fluctuation other than the expected Poisson vari- 
ation. Such a burst offers the best chance for a unified understanding of the gamma-ray 
burst and afterglow emissions. We provide a detailed modeling of the observed radia- 
tion from GRB 021211 both during the burst and the afterglow phase. The consistency 
between early optical emission (prior to 11 minutes), which presumably comes from re- 
verse shock heating of the ejecta, and late afterglow emission from forward shock (later 
than 11 minutes) requires the energy density in the magnetic field in the ejecta, ex- 
pressed as fraction of the equipartition value or e^, to be larger than the forward shock 
at 11 minutes by a factor of about 10 3 . We find that the only consistent model for the 
gamma-ray emission in GRB 021211 is the synchrotron radiation in the forward shock; 
to explain the peak flux during the GRB requires e# in forward shock at deceleration 
to be larger than the value at 11 minutes by a factor of about 10 2 . 

These results suggest that the magnetic field in the reverse shock and early forward 
shock is most likely frozen-in-field from the explosion, and therefore a large fraction 
of the energy in the explosion was initially stored in magnetic field. We can rule out 
the possibility that the ejecta from the burst for GRB 021211 contained more than 10 
electron-positron pairs per proton. 



Subject headings: gamma-rays: bursts, theory, methods: analytical - radiation mecha- 
nisms: non-thermal - shock waves 



1. Introduction 

Considerable progress has been made in the last few years toward an understanding of the 
nature of the enigmatic Gamma-Ray Bursts. Much of this progress has resulted from the observa- 
tion and analysis of afterglow emission - the radiation we receive after the high energy gamma-ray 
photons ceases - which is firmly established to be synchrotron radiation from a relativistic, external 
shock (e.g. Wijers, Rees, and Meszaros 1997). The nature of the explosion and the process that 
generates gamma-ray photons continue to be debated, although it is widely accepted that these 
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explosions involve a stellar mass object. The detection of narrow emission lines (e.g. Greiner et al. 
2003) and the emergence of a spectrum similar to that of SN 1998bw (Matheson et al. 2003) in 
the afterglow of GRB 030329 indicate that at least some GRBs are produced when a massive star 
undergoes collapse at the end of its nuclear burning life. 

Further progress toward understanding the GRB explosion requires afterglow observations 
at times closer to the burst and a simultaneous modeling of both the afterglow and gamma-ray 
emissions. In this way, one can explore distance scales of ~ 10 16 cm from the center of explosion, 
i.e. an order of magnitude smaller than that probed by afterglow emissions at half a day or later. 
Such a treatment is more likely to succeed in those cases where the prompt (burst) emission arises 
from the same region as the delayed (afterglow) emission, i.e. from an external shock. The simple 
FRED-like (fast rise, exponential decay) light-curves seen in about 10% of bursts represents the 
type expected from an external shock (Meszaros & Rees 1993), while short variability timescale 
bursts with complicated light curves are usually attributed to internal shocks in an unsteady outflow 
(Rees Sz Meszaros 1994; see Piran 1999 & Meszaros 2002 for recent reviews). 

This paper is an attempt to explain with the same process - synchrotron and inverse Compton 
emission from an external shock - the burst and afterglow emission of GRB 021211 detected by 
the HETE II (Crew et al. 2003), a burst which had a simple, FRED-like morphology and whose 
afterglow has been followed starting from 60 seconds until 10 days after the burst. In §2 we 
summarize the observations of GRB/afterglow 021211. In §3— §5 we present the formalism for 
calculating the synchrotron and inverse Compton emissions from both the forward and reverse 
shocks, and in §6 and §6.2 we assess the ability of the synchrotron, self-Compton model with a 
uniform and an r~ 2 stratified medium to accommodate the properties of the GRB 021211 and its 
afterglow. 



2. Summary of Observations for GRB 021211 

At a redshift z = 1.0, GRB 021211 had a duration of <; 2.3s in the 30-400 kev energy band and 
a fluence of ~ 2 x 10~ 6 erg (Crew et al. 2003); the duration in the 10-25 kev band was ~ 4s. If both 
the burst and the afterglow for GRB 021211 arise from some combination of reverse and forward 
external shocks, then the deceleration time is close to the time when the GRB light-curve peaks 
in 30-400 kev band, i.e. about 2 seconds. The average flux during the first 2.3s was 4 mJy in the 
7-30 kev band, 3 mJy in the 50-100 kev band, and 0.5 mJy in at 100-300 kev, while the peak of the 
vF v spectrum was at 47+y kev. The isotropic equivalent of the energy released in the 10-400 kev 
emission is ~ 10 52 erg. 

At 90s after the burst, the R-band magnitude of the afterglow 021211 was 14.06 (Wozniak et 
al. 2002), corresponding to a flux of 7.2 mJy. The optical flux decayed as t- l - 82 ± - 02 for the first 
10 minutes, after which it flattened to a i- - 82 ^ - 11 fall-off (Li et al. 2003), reaching magnitude 25 
at 7 days (Fruchter et al. 2002). The steeper decay seen during the first 10.8 minutes suggests 
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that the optical emission is dominated by the reverse shock energizing the GRB ejecta, while the 
shallower, later time decay is attributed to the forward shock that sweeps-up the ambient medium. 
The R-band flux at If minutes, when the two contributions are equal, is 0.39 mJy, therefore the 
forward shock optical flux at If minutes is 0.f9 mJy. Fox et al. (2003) report a 3-a upper limit of 
flO /j,Jy on the radio (8.5 GHz) flux at O.f days, and an upper limit of 35 fiJy during 9-25 days. 

Finally, Milagro has reported an upper limit of 4 x f0~ 6 erg cm -2 on the 0.2-20 Tev fluence 
over the burst duration reported by the HETE WXM (McEnery et al. 2002). 



Consider an explosion where the isotropic equivalent of energy release is £ and the initial 
Lorentz factor (LF) of cold baryonic material carrying this energy is Fq. Before the ejecta are 
significantly decelerated, the thermal LF j p j of the protons in the forward shock (FS), equal to 
the the bulk LF of the swept-up medium, is (e.g. Piran 1999) 



where n e j is the comoving particle density of the ejecta, and n = {A/m p )r~ s is the radial profile 
of the external particle density (s = for a homogeneous medium, s = 2 for a pre-ejected wind). 
The above result holds for n e j/n Tq, otherwise j p j « To- 

Taking into account that the laboratory frame energy per FS-heated proton is IVy Pj / = 7p/> 
the deceleration radius R d at which the energy of the swept-up medium is half the initial energy of 
the ejecta is 



As long as the distribution of LF of the ejecta is not too narrow, and the duration of the central 
explosion is less than R/(cTq), the comoving width of the material ejected in the explosion is 
proportional to R/Tq. Let us parametrize the comoving thickness of the ejecta as t]R/Tq. Therefore, 
the comoving density of the ejecta is 



3. Shock dynamics & deceleration time 



lpJ = r d * (r /2) 



1/2 (n e » 1/4 , 
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which substituted in equation (2) leads to 



(3 - s) 2 7]£ 
4ttAc 2 T 2 



i 

3-s 



R d = c 



(4) 



This result is identical to that obtained for the Blandford-McKee self-similar solution extrapolated 
back to R d if we set rj = 17/18 for s = and r] = 9/2 for s = 2. 
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From equations (3) and (4), the comoving density of the ejecta at Rd is 



-^-{Rd) — 7^ To^i ( 5 ) 

n (3 — s) 2 ?7 2 



for rj ^ 1. Substituting this into equation (1), the LF of the shocked ISM at Rd is 

r r Q 

d [2(3 

From equations (4) and (6), the observer-frame deceleration timescale is 



(6) 



t d -{l + Z)j v -^2 ~ rV 2 
ci d C 1 



(3 - sf- s rt*- s £ 



4ttAc 2 T 2 



i 



3-s 

(7) 



where f v is a correction factor that takes into account the difference between the arrival time 
(1 + z)Rd/{2cT 2 d ) of photons emitted from the contact discontinuity 1 and that from where most of 
the GRB emission arises 2 . The lab frame speed of the reverse shock relative to the back-end of 
the shell, shown in figure 1, is V rSt i a b/c f» 1.4(3 — s)rj /Fq. Thus, the time it takes for the RS to 
cross the shell (in lab frame) is, riRd/ (FQV rSt i a b) « [1-4(3 — s)]~ 1 Rd/c. And so the RS crossing time 
is same as the deceleration time to within a factor of order unity. 



For t > td the LF decreases as 

r(t) = ( 



<V>-r j{ lY"\ (8) 



4. Forward Shock 

The comoving density behind the forward shock (FS) is, p = 4poF, and the thermal energy 
density is u = ApQC 2 T 2 ; where po = Ar~ 2 is the density of the medium just ahead of the shock, 
and 7 is the bulk LF of shocked fluid given by equation (8). A fraction e e of the thermal energy of 
the shock- heated circumburst medium is taken up by electrons. Electrons with thermal LF greater 
than 7j are assumed to have a powerlaw distribution with index p, i.e. dN e /d'y oc for 7 > 7$, 
where 7$ = e' e (m p /m e )jp, is the minimum thermal LF of electrons; e' e = [{p — 2)/(p— l)]e e for p > 2 
and 7 P is the proton thermal LF. The energy density in magnetic field is assumed to be es f u, and 
therefore the magnetic field is B = 4Tc(27reB f Ar~ 2 ) 1 / 2 . 



lr The usual factor 2 in the denominator, corresponding to photons moving along the direction observer-center of 
explosion, is compensated by that most emission arises from the gas moving at an angle F^ 1 relative to that direction. 

2 For instance, if the burst is FS synchrotron emission from higher energy electrons in a fast cooling regime, then 
the 7-ray emission arises from the shocked gas immediately behind the FS and f v = 1/4. At the other extreme, 
when the burst arises from fast cooling electrons located immediately behind the reverse shock, it can be shown that 
/„ ~ (1/2)07/(3 -s) for V » 1. 
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The FS synchrotron injection frequency, i^/ and the flux at the peak of the F v spectrum, are 
qB^T Aqe' e 2 ei f mlA^T A f s r4 f 4 _ ,) rtjl -»/2 / t \ ~ 3 /2 



2irm e c(l + z) 



2tt m z A\ + 
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(10) 



d'l m p m e c(3 — s)d' L 2 

where q & m e are electron charge and mass, m p is proton mass, and d' L = Dl/\/1 + z, Dl being 
the luminosity distance, N e = AR ? '~ S / (3 — s) is the number of electrons per unit solid angle behind 
the shock, and e' u , = 3 1 ^ 2 q 3 B/m e c 2 is power per unit frequency per electron, in comoving frame, 
at the peak of the synchrotron spectrum. 

The synchrotron injection frequency for the cases of s = & 2 are written out explicitly for 
ease of application later on 



"if(t) = e' e 2 el f4 (t/t d ) 
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where A* = A/5xlO n g cm 1 , and an integer subscript n on a variable X, X n , means X/10 n . The 
flux at the peak of the synchrotron spectrum for s = & 2 is 



*/-4 



^2 



'5.5 x 10- 7 ^ /2 r« 2 *d mJy 



1.2 x 10 3 A^ /2 (£/i d ) 



s = 
1/2 mJy s = 2 



(12) 



In the derivation of the above equation we have set z = 1 which corresponds to the redshift of GRB 
021211. 

The FS synchrotron self-absorption frequency (vAf), obtained by equating the intensity at uaj 
to 2m e iifv\f, is given by 

' X " (6^) 1 /2^| r 4-3 S(4 ^ ) (2-3 S )/2 ( ^ 



; ,2 ^/ 

•A/ 



(13) 



^/ / 8(3 - s)m e m24(l + z )(6-3s)/2 c 3 S /2 \ td/ 

where a depends on the relative location of v A f with respect with v^f and the cooling frequency 
v c f\ for v c f > v if > v A f, a = 1/3, and a = p/2 if v if < v Af < v cf . 



4.1. Application to GRB 021211 late time optical observation 



We make use of the optical R-band flux at late time, t > llmin, to provide constraints on the 
density, e^/ & for s = and s = 2 cases separately. 
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4-1.1. Parameters for a uniform ISM model 
Using equation (7) the bulk LF at deceleration time, for s = 0, is found to be 



r 



<i 2 
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(14) 



where no is density of the uniform ISM, f v = max{l,2y^/3} (see footnote 2) and an integer 
subscript on a variable X n means X/10 n . 

The observed R-band flux at 11 min for GRB 021211 is 0.4 mJy. According to the fit presented 
in Li et al. (2003) the contributions from the reverse and forward shocks to the observed R-band 
flux are equal at this time. Therefore, the FS peak flux at 11 min is greater than 0.2 mJy; we will 
consider the peak flux to be 0.2 Af mJy, with Af > 1. It should be noted that for s = the peak 
flux is time independent. Substituting this into equation (12), and making use of equation (14) to 
eliminate we find 

4%*B,_ A = A ffv 3 - ( 15 ) 

The R-band lightcurve is observed to be monotonically declining from the earliest time (90s), and 
from 11 min to 10 days the decline is a simple powerlaw with index 0.82±0.11. Thus, the frequency 
of the peak of the spectrum at 11 min is expected to be less than the R-band frequency of 4.7xl0 14 
Hz or 1.95 ev. Let us assume that the peak frequency at 11 min is a factor A v smaller than the 
R-band frequency. Substituting this into equation (11) and making use of equation (14) we find 

^t/'e 2 /? 2 = 8xrt; 1 or e Bf (t = llmin) = 6AxW' 9 A^£ 5 - 2 1 e' e - i f- 3 . (16) 

The synchrotron peak frequency as a function of time is given by 

v if = 37.3t~ 3/2 A- 1 kev, (17) 

Combining equations (14), (15) and (16) we find 

T% « 2.2 (A v A f )-V*fW£ll 2 tt ^ ■ (18) 
Substituting this back into equation (14) we obtain 

n = 1.5 x 10 4 (A v AffS^f- 3 e' e A cm" 3 . (19) 



Note that Af and A v are related by Af = A^ 1 ^ 2 , if the cooling frequency (v c ) is above the R-band 

if = A V J 2 if v c is below the R-b; 
find n = 61 cm" 3 , T d = 250, and e B = 8 x 10~ 6 . 



at 11 min, and Af = A V J 2 if v c is below the R-band. For £52 = 1, e e = 0.5, A v = 4 & p = 2.2 we 
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4-1-2. Parameters for s = 2 model from late time Optical data 
The deceleration time for s = 2 is given by (see eq. 7) 

t d = 0.08(1 +z)^t- s, (20) 

where A* = ^4/5xlO n g cm -1 h f v = max(l, y/2rj/3). Since the peak flux at 11 min is 0.2 Af mJy, 
we find using equation (12), the FS peak flux at an earlier time to be 

f UpJ = 5AA f t- 1 / 2 mJy, (21) 

and 

ATeXji 2 V\= 4.5x10^ (22) 

Taking the synchrotron peak frequency at 11 min to be 4.7xl0 14 vl~ 1 Hz, and substituting this into 
equation (11) we obtain for s = 2 

Al^lTl/J 2 = 5 x 10- 3 47\ (23) 

and the time evolution of v^f is same as in equation (17). 
From equations (22) and (23) we obtain 

A, = 0.9< 2 (4/^/), (24) 
which when substituted into equation (20) gives the LF at deceleration 

r d2 = o.T^V^V-v^/r 174 - (25) 

Combining equations (23), (24) and (25) we obtain 

e Bf (t = 11mm) = 1.7 x 10 -8 A' 2 f^S^e'^ 4 . (26) 

We see from these equations that for s = 2, A v can be as large as 10-20, and yet give acceptable 
values for various parameters. 



5. Reverse Shock 

The emission from reverse shock (RS) in gamma-ray bursts is discussed by a number of authors 
e.g. Panaitescu & Meszaros (1998), Sari & Piran (1999), Kobayashi (2000), Piran (2000). A 
particularly important parameter that determines the behavior of RS is the thickness of the shell 
of material or ejecta that carries the relativistic energy of the explosion. We have parametrized the 
ejecta thickness as r}RjT\ in lab frame; for a shell whose thickness is dominated by expansion at 
deceleration radius we expect i] ~ 1, otherwise the thickness is determined by the duration of the 
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central engine, and 77 could be much larger than unity at ta- We calculate RS emission for a range 
of 7] between 0.5 and 10. Fortunately the main conclusions of this work for GRB 021211 remain 
unchanged even for a larger range of rj. 

At the deceleration radius R d , the ratio of the thermal energy of protons in the reverse shock 
(RS) region, j Pjr , to that in the FS, j p j, is (see fig. 1) 

Tfer ^ J_ fn^Y°- 7 „ [(3~^] L4 _ [(3 ~ s) V ] - 9 (2?) 



ipj 4r o Vrgny 4r v^r^ 

The first part in the above equation is valid only for 0.01 ^ n e jjY\n ^ 100; for n e jjY\n <C 0.01, it 
can be shown that 7p,r/7p,/ — i n ej / n)~ ' 5 ■ In deriving the second part of this equation we made 
use of (5) for the density of the ejecta at Rd — n e j(Rd)- It should be noted that the thermal energy 
per proton in RS is ~ m p c 2 rf '- 9 /32 1 / 2 , and so protons are not heated to a relativistic temperature 
in the reverse shock. 

The pressure continuity across the contact discontinuity surface, which separates forward and 
reverse shocks, implies that the magnetic field strength in RS and FS are equal, provided that sb 
is the same behind both shocks. However, one might expect e# in the RS (esr) to be different from 
the value in FS (ee/). Then the synchrotron peak frequency in the RS is 



This can be written out explicitly as follows 



* rldj " (1287T)V2 ( l + z)m 3 ' ^ 



or 



.*(«=4^4 3 (^) l ^ l ■ 8 x( l ■ lxloI l 6 x /4Hz s =° (3o> 

U.lxlO 15 ^ 1 Hz s = 2 

In deriving this last equation we made use of equations (18), (19), (24) & (25) for ISM density and 
LF at deceleration. 

Since the FS and RS region are moving at same LF, at deceleration, the RS peak synchrotron 
flux is equal to the FS peak flux times the ratio of number of electrons in the ejecta to the swept-up 
electrons in the surrounding medium up to R s ; this ratio is equal to Tq/(3 — s)i] = [2/(3 — s)^] 1 /^^. 
Thus, the RS peak flux is 



fu p ,r(td) - fu p ,f(td) 



(3 - s)ij 



1/2 



r d , (3i) 
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or 



(3e Br A)^q 3 £{l + z) 
Using equations (18), (19), (24) & (25) this equation reduces to 



fv p ,r(td) — . i- (32) 

m e m p c3d>lr d -\4ct d ) a [vr(3 - s) V ] a 



/^) = 4\4 3/2 (W^ mJy 8 = (33) 

U.8xl0 4 t d d/ mJy s = 2 



The RS synchrotron self-absorption frequency is 

o fn f%r_ g 3 ^ [6^(1 + z) s ] l/2 (u) 

v AA t d) — - , r , ; s +4„ , > (34 j 



>"W 8vrm e m2 c 3 e / e [(3 - s)r/] L4 (4ct d ) £ ^r* +2 ' 

where a = 1/3 if fj r > VAr, and p/2 otherwise. Using equations (18), (19), (24) & (25) this can be 
rewritten as 

" 4 U.lxlO 15 ^ 1 Hz s = 2 



5.1. Time dependence of radiation from reverse shock 

The time dependence for Vi r & /i/ p ,r is determined by the evolution of magnetic field and 
electron thermal energy in the reverse shock. Electrons in the ejecta cease to be heated after 
the passage of the RS, and their energy decreases with time as a result of adiabatic expansion. If 
electrons continue to exchange energy with protons, and the fraction of thermal energy in electrons, 
e e , is time independent, then electron thermal LF decreases as 

7 e oc {R 2 5RY^ oc t'^^T, (36) 

where SR is the comoving shell thickness which is a weak function of time for sub- or mildly- 
relativistic RS, and R, the radius of the ejecta, increases with time as i 1 /*- 4 "^; e e = 1 if electrons 
and protons are decoupled. 

The magnetic field, frozen in the ejecta, decreases as B' oc (RSR)^ 1 oc t _1 /( 4_s ) if the field is 
transverse; a longitudinal field decreases as t~ 2 /( 4 ~ s ), therefore any non-zero transverse field will 
become the dominant component at large distances. 

The synchrotron injection frequency and the peak flux decrease as 

v ir oc B'j^T oc t~ 31 6(4-4 Efi ? and f Up r K N e B'Y oc N e t~^^ , (37) 
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where N e is the total number of "radiating" electrons in the ejecta. 

The cooling frequency decreases at the same rate as z/j r ; the decline is faster if radiative losses 
dominate over adiabatic losses. The number of electrons radiating in an observer band might have 
a non-trivial time independent if the magnetic field is not constant across the ejecta - electrons in a 
region of higher magnetic field will lose energy at a higher rate and their radiation drops below the 
observed band sooner than electrons in lower magnetic field region. This together with uncertainty 
with the evolution of 7 e - which depends upon coupling between electrons & protons - and the 
unknown energy density & LF structure of the ejecta, makes it difficult to calculate with confidence 
the power-law decay index for flux from reverse shock. 

In order to fit the data for GRB 021211 what we do instead is to work backwards from the 
observed lightcurve slope and determine the decay of synchrotron frequency - which depends on 
both B' and j e and so its time dependence is more uncertain than the peak flux which depends 
on B' alone - that we need to calculate the cooling frequency at 11 minutes after the explosion to 
make sure that it is above the optical R-band. 

Let us consider the time dependence for injection frequency and peak flux to be t~ a " & t~ a f 
respectively. The flux above the synchrotron peak decays as i~[ 2a /+(p _1 )°' y ]/ 2 (jp is electron energy 
power-law index). The observed decay for GRB 021211 was i~ L8 , which we use to determine a u ; 
we assume that a f is as given in equation (37), but allow for a small deviation when fitting the 
observed data. With a v and ctf thus determined, we find the time dependence of B' & 7 e 

B' oc t~ af+ ^ , 7e oc r^, (38) 
that we use to calculate absorption and cooling frequencies and flux as a function of time. 



5.2. Compton Parameter & Cooling Frequency 

The comoving frame timescale for an electron of energy m e c 2 j e to cool as a result of synchrotron 
and inverse Compton emission is 

, = 67rm e c = 3m e R s d 

c a T B'^ e (l + Y) 16a T ce Ble T 2 d A(l + YY { > 

where Y is the Compton Y parameter, prime denotes comoving quantity, and the cooling is con- 
sidered at the deceleration time. At deceleration, when t' c = 4*^^/(1 + 2), the electron cooling LF, 
7c defined by the equality of the radiative and dynamical timescales, is 

= 3vrm e c(l + z) = 3m e r^~ 3 

lc{d) 2a T B'H d T d (l + Y) 16a T e B A(l + zy-HAct d y- s (l + Y) { ) 
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Substituting for A and from equations (18), (19), (24) & (25) we obtain 

£ U 4 /513 s = 

lc(td) = — — x S (41) 

e B _ 4 e'J(A,A f )i(l + Y) \ 22t f s = 2 

The cooling frequency v c , defined as the synchrotron frequency for electrons with LF j c , is 

qB'jfrd _ , 1 x 5 c 2 (4ct d r 2 )^ 
" M ~ 2vrm e c(l + z) " ^ X 10 {A £B _J3/2 (1 + y )2(1 + z) fr-2),2 Hz ' ( 42 ) 

which can be rewritten by substituting for A and 

t . fa) - 4 /2 f2xl0 15 (l + Z )^ 2 Hz s = 

e l 2 4 6'4(A^ / )2(l + y)2 l 9 xl0 13 (l + z)- 2 4 /2 Hz s = 2 



The cooling frequencies in the reverse and forward shock regions are calculated from this 
equation by substituting appropriate values for e# and Y corresponding to each region. The 
Compton Y parameter is calculated below. 

The electron column density in the ejecta, at the deceleration radius, assuming that the ejecta 
consists only of protons and electrons, i.e. there are no pairs, is 

N = £ = : a + s) 2 * (44 ) 

e ' r 47rR 2 d T m p c 2 64^^2(3 - s)r] m p cH 2 d T d ' 1 ; 

The optical depth of the ejecta to Thomson scattering is 

(1 + z) 2 E 

T r = a T N er = 1.7xlO~ 3 f. -. (45) 

The Compton parameter Y = t~/ 2 , where 7I, the mean squared electron LF, for va < v % < v c 
Sz 2 < p < 3 is 

and m e c 2 7j is the minimum thermal energy of shock heated electrons. The Compton parameter in 
the particular case of va < v c < V{ is given by 

Y,(£ . (47) 
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Substituting (41) into (46) and making use of equations (27) & (45) we find the Compton Y 
in the reverse shock region 

Y(l + Yf-v = ^ — - 

(p- 2) (3- p)m{T 32— 

f 2.8 x 10 4 191"% 8 s = 



^6.4 x 10 2 22 _P ^ 4 s = 2 



d (48) 



The Compton Y substituted back into equation (43) yields the cooling frequency in RS. A 
similar calculation gives u c in the forward shock. 

When ua > min{f c , i/j}, the synchrotron photon flux that is scattered by an electron is dimin- 
ished by self-absorption and va and v c have to be determined by solving a set of coupled equations, 
as described in Panaitescu & Meszaros (2000). Some of the cases considered for GRB 021211r 
fall in this more complicated regime, and all of the numerical results presented in this paper are 
obtained by determining va and v c numerically, in a self-consistent manner. 



6. A unified modeling for 7-ray and afterglow data 

We apply the results of the last two sections to a systematic analysis of 7-ray, optical and radio 
observations for GRB 021211 and determine models that are consistent with all data. We discuss 
the cases of a uniform density ISM (s = 0), and a medium carved out by the progenitor's wind 
(s = 2) in two separate subsection. 



6.1. A Uniform Density Circumburst Medium (s=0) 

In the next subsection we discuss the early optical & radio emissions from the RS. We take 
up the question of what could have produced the 7-ray emission in §6.1.2 & radio upper limit in 
56.1.3. 



6.1.1. Optical and radio emissions from reverse shock 

The RS synchrotron injection frequency is a factor of (7 p ,//7 P ,r) 2 smaller compared with the 
peak frequency of the FS emission as long as 63 is the same in reverse and forward shocks. The syn- 
chrotron injection frequency in the RS, otherwise, is expected to be about 2xlO~ 2 ri 1 - 8 A~ 1 (eBr/£Bf) 1 ^ 2 
ev (see eqs. 30 & 16). This suggests that the RS flux in the optical band decreases with time for 
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t > td- The extrapolation of the observed flux of 7.2 mJy at 90s, with a powerlaw decline of i~ L8 , 
gives a R-band flux at 5s of 1.3 Jy or 8.5 mag. So, contrary to claims, GRB 021211 was as bright 
as 990123 close to the deceleration time. 

The injection frequency in RS for s = declines with observer time approximately as and 
the cooling frequency too declines as t~ l or faster. After the passage of the reverse shock, which 
takes place on the deceleration time scale of a few seconds for GRB021211, electrons are no longer 
accelerated and there is no emission from RS at a frequency greater than the cooling frequency 
(v cr ). The R-band flux from GRB 021211 is observed to decline as i -1 - 8 for 11 minutes and then 
the decline slows down to i -0-8 . This suggests that the RS emission lasts for at least 11 minutes in 
the R-band, and therefore v cr at deceleration should be ^ 10 17 Hz. 

The inverse Compton parameter and the cooling frequency are determined from equations (48) 
and (43). For p = 2.5 these quantities at the deceleration time are 

Y(l + Y) - 5 = 2 x 10 3 ef 75 €gJ rp M Ap£^t~^ , (49) 



and 



,1/2 



■ cr = 1.5xl0 11 b - ^ Hz, (50) 

/7.7.0.83 „1.1 ' ^ > 

e e B r A V A v l d 



S, 

so long as Y ^ 1. For Y <C 1 the cooling frequency is obtained by setting Y = in equation (43). 

The requirement that v C r{td) > 10 17 Hz - in order to have non-zero flux in the R-band from 
RS for ~ llmin - provides an upper limit on e' e given below 

175 fV 15A 

Substituting this into (30) & (33) we find the injection frequency and peak flux from the RS 

u ir (t d ) < 6 x 10 13 e^/ 5 f X^V^V' 9 Hz, (52) 

and 

/, P ,,(t,)^2194 / r 3 _ 4 ^ 2 85 ^ /20 r / -°-V- 45 mJy. (53) 
The flux in the R-band {i>r = 4.95 x 10 14 Hz) at deceleration time is 

U^^me^fE^/^hd 11 mJy, (54) 

whereas the flux in R-band at 90s is fu R (td){t/td)~ 1 ' 8 , if Vi r {td) < ^R, otherwise the flux is given by 

f VR (t = 90s) = U p Atd){Ws/t R y l -\t R /t d )- a f, (55) 

where t R is the larger of t d & the time when Vi r {t) = 4.95 x 10 14 Hz (R-band frequency), and a/ ~ 1 
is the power-law decay index for the peak-flux i.e. fv p ,r(t) oc t~ a f . The above equations for R-band 
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flux are applicable when the synchrotron-self-absorption frequency is less than vr which is indeed 
the case as (see eq. 35). 

Since v ir (t d ) < 4.95 x 10 14 Hz, the R-band flux at 90s is found to be 0.014 e^f^A^^r] - 3 ^ 7 
mJy. A R-band flux of 7.2 mJy, as observed for 021211, requires e_e r ~ 10~ 2 , £52 ~ 20 & A v ~ r\ 
(td ~ 3s). Figure 2 shows the allowed parameter space which satisfies the R-band flux at t > 90s. 
Note that e^/e^ ~ 10 3 , for the allowed parameter space. This result is consistent with the 
analytical calculation presented above. 

Fox et al. (2003) reported that the flux at 8.5GHz, 0.1 day after the burst, was less than 
110/t/Jy. This frequency is above the self-absorption frequency, and the flux from RS is a few 
times larger than this upper limit for the parameter space in fig. 2, for ~ 3s. Diffractive 
interstellar scintillation can decrease the flux at 8.5GHz at this early time by a factor of a few 
thereby providing consistency with the reported flux upper limit. We note that the radio flux 
would exceed the observational limit by almost an order of magnitude if the time for RS crossing 
were taken to be 30s, almost independent of the details of RS model, which suggests that the shock 
crossing time is approximately equal to the burst duration of 2-4 s. 



6.1.2. ~f-ray emission during the GRB 

The injection frequency in FS at the deceleration time (td = 2 s) is 13.2 A^ 1 kev (see eq. 17), 
and the peak flux is 0.2^4/ mJy. To explain the early optical afterglow requires A u to be of order a 
few or larger (see fig. 2 & the discussed in the last subsection), and therefore the injection frequency 
and the peak flux in FS at t^ are ~ 5 kev & 0.4 mJy respectively. The observed values for the peak 
7-ray flux during the GRB is about 4 mJy, and the peak of vf v is at ~ 50 kev (Crew et al. 2003). 
Thus, the observed peak flux during the GRB is about an order of magnitude larger than predicted 
by the extrapolation of the optical data at 11 min. And the observed peak frequency, depending 
on the value of the cooling frequency in forward shock, is also about an order of magnitude larger 
than the synchrotron peak frequency. 

We consider whether synchrotron-self-Compton process in the reverse or the forward shock 
might explain the gamma-ray emission. The peak of vf v for inverse Compton scattered synchrotron 
photons occurs at a frequency of max{7 2 , 7 2 } x max{z/j, u c } (see Panaitescu & Meszaros, 2000). For 
the reverse shock of GRB021211, v c ~ lkev, and 7 C ~ 3 x 10 2 (see §6.1.1), and therefore, the IC 
peak frequency is at ~ 10 2 Mev, or three order of magnitude above the observed peak; the IC flux 
at 50 kev is about 0.1 mJy. For the forward shock, V{ ~ 5kev h 7« ~ 10 4 , and thus the IC spectrum 
peaks at a energy ^ 0.5Tev; the flux at this energy is smaller than the upper limit provided by 
Milagro (McEnery et al. 2002). 

Having eliminated synchrotron-self-Compton process as an explanation for the 7-ray emission 
from GRB 02121, we turn to synchrotron emission from the reverse or the forward shock as a 
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possible mechanism to account for the observations 3 . 

The synchrotron emission from RS can have v c ~ 50kev provided that we consider a small 
value for e' e ~ 0.04. The flux at 50 kev can be calculated directly from the observed optical R-band 
flux and is estimated to be about 4 mJy - consistent with the observed 7-ray flux. This would 
have been a very economical and elegant explanation for all the observations for 021211 from 7-ray 
to radio frequencies. However, this possibility is, unfortunately, ruled out by the observed spectral 
slope of 0.4 (/„ oc u 0A ) below the 50 kev peak (Crew et al. 2003), whereas the RS synchrotron 
model predicts a spectral power-law index of — (p — l)/2 ~ —0.5. 

As we discussed earlier, synchrotron emission from the forward shock cannot account for the 
gamma-ray observations as long as we take e# in FS at deceleration to be same as it is at 11 
minutes. Having ruled out all possibilities for producing 7-rays in a standard external shock (for 
s = 0), we now relax the assumption of time independent e# in the FS. Our goal is to find a solution 
where v c ~ i/j ~ 50kev, and the peak flux is 4 mJy; note that v c ~ i/j is required by the observed 
low energy spectral index 4 . 

For a time independent €b in the forward shock Vif(td) ~ 37.3t d 3 ^ 2 A^ 1 Kev (see eq. 17), and 
the peak flux is 0.2 A f mJy. Therefore, the flux at 50 kev is ~ 0.2 Af(u if /50Y P ~^ /2 ~ 0.16 tT 9/8 mJy; 
note that Aj = A^ 1 ^ 2 . This flux is too small by a factor of about 10 2 to satisfy the 7-ray 
observation for 021211. The only way out of this difficulty is to consider e# to be larger at the 
deceleration time by a factor of ~ 100 than its value at 11 minutes. Since the flux at v \ < v < v c 
scales as e^ -1-1 ^ 4 , and the synchrotron injection frequency vi oc e]/ 2 , a larger e# by a factor 10 2 
will increase the flux at 50 kev by two orders of magnitude & increase v^j by a factor 10, thereby 
simultaneously satisfying the observed 7-ray flux and the peak frequency requirements. The rest 
of this section is devoting to the calculation of cooling frequency in FS at the deceleration time to 
make sure that a larger value of cb at td is compatible with the requirement of u c f ~ 50kev. 

The optical depth to Thomson scattering in forward shock at the deceleration time is 



°Tn G R d _ 0f . ^_ & e' e \A f A v f/H 



1/4 
d 



= j. » a = g g 1Q - b g v j v, d , 5g ^ 

f 3m p (i + ,)4 /2 

where = 4ct,iT d / (I + z) is the deceleration radius, and the last equality was obtained by making 
use of equations (18) & (19) to eliminate (LF at deceleration) & no (ISM density). The Compton 
parameter Y = Tfy% with 7 2 given by equation (46) when j c > 7,; 7, = (m p /m e )e' e Td in the forward 
shock at td- Combining these equations we find 

-y— r/=4 . 5xl0 -a <p-wt\ ^y-'4^(w^r, (57) 



if (3 -p)(p- 2) ' (3-p)(j,-2) \mj e \n 



3 We feel that a simple single peaked lightcurve for 021211 should not require internal shocks to produce 7-ray 
emission, which were invoked to explain multi-peaked and highly fluctuating GRBs. 

4 For any other ordering of v c & z/» the spectral index below the peak of vf v will be close to -0.5. 
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Substituting for 7 C from equation (41) we obtain 

B f-4 V d ) 5 



"52 

For p = 2.5 and 7>1 the Compton parameter is 

F = 2xl0 3 e B f 4 4t(^)h- 1/6 . (59) 
Substituting this into equation (43) we find the cooling frequency in forward shock for z = 1: 



^c/ = 7xl0 9 - 14 - Hz. (60) 

When i/ c/ < v { & y » 1, Y « (e e /e jB/ ) 1/2 , and 

•"> = 2xWll ^\ + M" Hz ' <61> 

e B/_ 4 Cee e ^ c (2 

For y <C 1 the cooling frequency in FS can be obtained by setting Y = in equation (43). 

The peak of the spectrum for 021211 is at v c f ~ Vif = 50kev = 1.2 x 10 19 Hz. Substituting 
this into equation (60) gives 

e' e A¥ej f Jt d ) = 0.055 £# tj*. (62) 

Combining this with equation (16) - under the assumption that e' e is time independent - and 
making use of the requirement that €B f (td)/eB f (t = 11mm) ~ 10 2 discussed earlier, we find 



z^USM =A3A?tf. (63) 

This relation can be satisfied if we consider, for instance, £52 = 20, €B f {td) ~ 10 -2 k, A v ~ 1, and 
provides a self consistent solution that accounts for 7-ray and optical radiations for GRB 021211. 

Figure 2 shows the parameter space allowed by the 7-ray flux for 021211 originating in the 
forward shock. Note that there is a range of parameters for which the early and late optical, and 
7-ray observations can be simultaneously explained. The general requirement, however, is a large 
magnetic field in the ejecta, somewhat smaller field in the forward shock at deceleration, and a much 
smaller es f at 11 min when the forward shock emission starts to become a dominant contributor 
to the optical flux. 
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6.1.3. Radio flux limit at 10 days 

The flux for GRB 021211 at 8.5 GHz, 10 days after the burst, is reported to be less than 
35 ujy (Fox et al. 2003). The synchrotron peak frequency, which decays as i~ 3 / 2 , independent of 
ISM density stratification, is 9.7 A^ 1 GHz at 10 days. If the radio band frequency were above the 
synchrotron self-absorption frequency, then the flux at 8.5 GHz should be 0.19 mJy independent of 
A u , which is a factor of 5.4 larger than the observed upper limit. 

We consider the possibility that the self-absorption frequency is larger than 8.5 GHz by a factor 
of about 2 thereby reducing the flux below the observational upper limit. The FS self-absorption 
frequency, calculated using equations (13), (16)-(19), is 

v A Mv A )V* = - ™*)^,A * A n o 4t Tt~ 1/4 Hz = 2.4 x 10^/ 2 £-VVVVV4 Hz. (64) 
(1 + z) 6 ' 2 e e 

In deriving the second equality we have set z = 1. Thus, va ~ 1-2 x 10 8 Hz at 10 days for A v = 4, 
£52 = 10 & e' e = 0.03. We see that A v should be ~ 100 in order that va ~ 20 GHz, and the flux in 
8.5 GHz band at 10 days is below 35 ujy. However, for A v ~ 100, n ~ 2 x 10 6 cm" 3 & e B ~ 10~ 8 
(see eqs. 18 Sz 19); these values are in contradiction with the requirement that no ~ 10 -2 cm -3 & 
CBf(t = llmin) ~ 10~ 4 in order to produce 7-ray emission in external shock (see fig. 2). 

A decrease in the density with r, such as for s = 2 medium, can reconcile the radio flux limit 
at 10 days. However, it is very difficult to produce 7-ray radiation, as seen in 021211, for s = 2 
density stratification (see §6.2.2). 

A loss of explosion energy by a factor of about 10 between 11 min and 10 days would also reduce 
the radio flux to a value below the observational upper limit (the peak flux & synchrotron peak 
frequency are proportional to £ & £}l 2 respectively). This requirement is, however, inconsistent 
with small e# at late times (t > llmin), and the fact that u c > i/j at all times. 

Yet another way that the radio flux at 10 days can be reduced is by requiring the jet break 
time to be less than 10 days. The isotropic equivalent of energy in 021211, estimated from early 
optical data (§6.1.1), is 2xl0 53 erg, which suggests jet opening angle to be about 8° (energy in 
GRBs is found to be narrowly clustered around 10 51 erg - Panaitescu Sz Kumar 2002, Frail et al. 
2002, Piran et al. 2002). For no ~ 10~ 2 cm -3 (see fig. 2), we expect the jet break time to be about 
10 days - which is consistent with the lack of a clear break in optical light-curve - and this can 
reduce the radio flux by a factor of a few. Another factor of 2 decrease could come from a decrease 
in £b by a factor of 2 between 11 minutes k, 10 days 5 . These effects together could then reconcile 
the late time radio flux in the model considered here with the observational upper limit. 



5 The decline of optical lightcurve as t 0,82 together with the optical spectrum of v 9 limits the decline of ta 
between llmin and 10 days to be less than a factor of ~ 3. 
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6.1.4- Milagro limit 

Milagro reported an upper limit of 4 x 10 -6 erg cm -2 on the 0.2-20 Tev fluence for GRB 021211 
(McEnery et al. 2002). We discuss whether the solutions we have found for early optical and 7-ray 
radiations, shown in fig. 2, respect this limit. 

The peak of the inverse-Compton in reverse shock is at f c 7 2 ~ 10 2 Mev, and the fluence in 
Milagro band is estimated to be ~ 10" 8 erg cm~ 2 . The IC emission in the forward shock peaks at 
Vijf ~ lOTev, and the fluence in 0.2-20 Tev is ~ Yvif Ui t d ~ 5 x 10~ 7 erg cm -2 - smaller than the 
Milagro upper limit by an order of magnitude. 

6.1.5. A summary of results for a uniform- density medium 

A uniform density medium cannot simultaneously explain the R-band afterglow emission after 
11 minutes (FS emission), and before 11 minutes (presumably from the RS) unless the energy 
density in magnetic field in the RS is at least a few hundred times larger than the magnetic energy 
density in FS. Moreover, for s = and an external shock origin for GRB 021211, the 7-ray emission 
can be produced via the synchrotron process in the forward shock provided that the magnetic field 
parameter (e#) in FS at deceleration is larger by a factor ~ 10 2 compared with the value at 11 
minutes, i.e., it requires a time dependent es f during the first few minutes. To satisfy the upper 
limit on radio flux at 10 days seems to require a combination of jet break at about 10 days and a 
decline of by a factor of ~ 2 between 11 min and 10 days. 

6.2. Pre-Ejected Wind Circumburst Medium (s=2) 

We consider early optical emission from reverse shock in the next sub-section, 7-ray emission 
in §6.2.2, and radio data in §6.2.3. 

6.2.1. Optical emission from reverse shock for s=2 

The synchrotron injection frequency and the flux at the peak of the RS spectrum are given by 
equations (30) & (33). For p = 2.5, £ 52 = 10, A v = 4, A f = A^~ 1)t2 = 2.8, t d = 2s, 77 = 4, and 
e' e = 0.05, we find v ir = 1.3 x lO^e^Hz, and f Up>r = 7Ae^_Jy; note that T d = 234, A* = 0.03 & 
CB f (t = 11mm) = 1.7 x 10~ 5 for this choice of parameters. The resulting flux in the R-band at 90s 
is ^3 0.5mJy or a factor 15 smaller than the observed value if we take €B r = 1- For e# r ~ 10 2 the 
flux agrees with the optical observation, provided that absorption frequency is below the R-band 
and the cooling frequency is sufficiently high, v cr ~ 10 17 Hz, so that electrons in the RS continue 
to radiate in the R-band for ~ 11 minutes. We look into these requirements below, and determine 
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the set of parameters that satisfies optical observations. 

The RS synchrotron self-absorption frequency at deceleration time is given in equation (35). 
For the parameters considered above the self-absorption frequency is ~ 10 14 Hz, below the R-band, 
and does not affect the optical flux. 

The calculation of cooling frequency proceeds in the same manner as for the s = case 
considered in the previous section. The inverse Compton parameter and the cooling frequency are 
determined from (48) & (43). For p = 2.5 & s = 2 these quantities at the deceleration time are 

Y(l + Y)™ = 1.1 x 10 4 4 2 - 75 e - r 1 / 2 4 r / - 85 4 5 /32^-i/8^3/8 ) (65) 

and 

E 1/2 t, 

"cr = 9 ' 5 >< 10? ,7.7 e 0.83 ZTTTB Hz > ( 66 ) 

so long as Y ^ 1. For Y <C 1 the cooling frequency is given by setting Y = in equation (43). 

The requirement that v cr (td) > 10 17 Hz - in order to have non-zero flux in the R-band from 
RS for ~ llmin - provides an upper limit on e' e 

£ e~ 1/9 ,2/3 1/7 • 

Substituting this into (30) & (33) we find the injection frequency and peak flux in the RS 

u tr (t d ) < 1.5 x lO^ ^X-' V ^ - 6 Hz, (68) 

^ 660 4 3 _ 4 ^ 85 ^ 2 %- a7 i d -°- 55 mJy, (69) 
and the flux in R-band at deceleration time is 

f VR (t d )< 8.5 e^^A- 7 ^^ 1 - 1 mJy. (70) 

The flux in R-band at 90s is therefore about 2.6 x 10~~ 3 e%f £ 52 Av 7/16 T]°- 3 t d - 7 mJy. This is a factor 
2 smaller than the observed value of 7.2 mJy, even when we takee Br _ = 10 4 , £ 52 = 20, and A v = 1. 
So, to obtain the desired optical flux at 90s, in a pre-ejected wind circum-burst-medium model, 
requires very extreme, and perhaps unphysical, parameters. More accurate numerical calculations 
support this conclusion. 



6.2.2. 7 -ray emission for s=2 medium 



The arguments against synchrotron self-Compton in reverse or forward shock as an explanation 
for the 7-ray observations for GRB 021211, we suggested for s = 0, also apply to s = 2. So we 
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once again turn to synchrotron radiation in forward shock, the most likely mechanism for 021211, 
to explain the 7-ray observations. 

For a time independent en in the forward shock the synchrotron injection frequency at the 
1 3/2 

deceleration is 37.3 A~ t d kev - a factor ~ 4 smaller than the observed v pea k ~ an d the peak flux 
is 5A/tJ 1/2 mJy. The flux at 50 kev is ~ 5A f t d 1/2 (z/i//50) (p ~ 1)/2 ~ 4i~ 13/8 mJy, which is also a 
factor ~ 4 smaller compared with the observed flux. Removing these discrepancies requires €b s at 
td to be larger by an order magnitude than the value at 11 minutes. 

As discussed in last section we require v c f ~ v^j ~ 50 kev in order to be consistent with 021211 
7-radiation. We compute the cooling frequency below to determine whether this condition can be 
satisfied for s = 2. 

The equation for Compton parameter in FS for s = 2, when 7 C > 7$, is derived in the same 
way as the s = case considered in §6.1.2, and is found to be 



Y(l + Yf-v = 7 4 x 10-4 fr- 1 ) teV" ^(AfA^r?- 1 
{ ' " (3-p)(p-2) \mj s-W^.s/^p-s ^ i^' 

e B f 4 \ ZjZl d ) c 52 l d 



(71) 



For p = 2.5 and F>1 the Compton parameter is 



Y = 8.7 x 10 3 e B ; /3 4 e; 5 / 3 ^/ 12 t d 1/3 . (72) 
Substituting this into equation (43) we find the cooling frequency in forward shock for z = 1: 

cl/2,7/6 

"c/ = 1-5x10 ^ g 52 22 ; 344/3 Hz. (73) 



For i/ c / < i/j & y » 1, y « (ee/es^ 1 / 2 , and 



ci/2,1/2 



B f-4 e e 



Hz. (74) 



For y <C 1 the cooling frequency is obtained by setting Y = in equation (43). 

Since the peak of the 7-ray spectrum for 021211 is at v c f ~ v%f = 1.2 x 10 19 Hz, we obtain from 
equation (73) the following relation 

e' e A?ef f (t d ) = 8.4 x 10- 3 £ 5 f tf . (75) 

Using equations (26) and (75), and taking eB f (td)/^B f (t = 11mm) ~ 10 as discussed at the 
beginning of this subsection, we find 

e%(t d )£itf =33Ap. (76) 
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This equation can be satisfied if we take £52 = 30, 6B f (td) ~ 0.1, td ~ 2s & Aj, ~ 1. However, for 
these parameters, the density of the medium A» ~ 5 x 10~ 4 , which is three orders of magnitude 
smaller than the value for a typical Wolf-Rayet star wind. Thus we find that the s = 2 model has 
difficulty producing the early optical and 7-ray flux for 021211. 

6. 2. 3. Radio flux upper limits 

The FS peak flux for 021211, for s = 2 medium, declines as 574/i~ 1//2 mJy, and therefore the 
peak flux at 10 day is ~ 5.8/uJy. The expected flux at 8.5 GHz at 10 day after the explosion is thus 
less than 5.8/xJy, and entirely consistent with the upper limit of ~ 35/uJy (Fox et al. 2003). 

The peak frequency for the RS emission declines as i~ 17 / 12 and the peak flux declines as ~ t -1 . 
Therefore the peak frequency and flux at 0.1 day are 8 x 10 8 Hz and 1.0 mJy, respectively. The 
absorption frequency decreases as v\ oc t~ 5 / 6 and thus ^ ~ 1.1 x 10 11 Hz at 0.1 day. Therefore, 
the flux at 8.5 GHz at 0.1 day is expected to be about 6^Jy, which is well below the upper limit of 
Fox et al. (2003). 

6.2.4- A summary of results for s = 2 medium 

Pre-ejected wind circum-burst medium (s = 2), for 021211, very easily accommodates the flux 
upper limit in 8.46GHz band reported at 0.1 and 10 days. However, unlike, s = 0, it requires 
extreme parameters ~eB / ~l&A^~5x 10~ 4 - to explain the 7-ray & early optical radiations. 
We did not find any common set of parameters, at 11 min, to explain both the early afterglow 
and the 7-ray observations simultaneously; solutions to explain early optical emission require, as in 
s = 0, magnetic field energy density in RS to be larger than the FS by a factor of a few hundred. 

7. Discussion 

We have shown that for GRB 021211 the gamma-ray burst emission, the optical afterglow, 
and upper limit on radio flux, can be understood in terms of emission from forward and reverse 
shock regions. To explain the early optical data - prior to 11 min - as resulting from the reverse 
shock, requires a high energy density in magnetic field, es r ~ 10 _1 , which is about three orders 
of magnitude larger than the value we obtain for the forward shock region at 11 min from optical 
data. Zhang et al. (2003) have suggested that the bright optical flash in GRB 990123 also requires 
high magnetic field in reverse shock. 

We find that the gamma-ray emission for GRB 021211 can be explained as synchrotron radi- 
ation in the forward shock. The value for eb in forward shock at deceleration required to explain 
the GRB fluence is about 10~ 2 which is larger by about two orders of magnitude compared with 
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the value at 11 minutes. This suggests that e# in the forward shock is time dependent at early 
times. The combination of high e# in the ejecta, and high but somewhat smaller in forward 
shock at deceleration time, suggests that the magnetic field in the RS is perhaps the frozen-in field 
of the highly magnetized ejecta from the burst as suggested by eg., Usov (1992), Meszaros & Rees 
(1997), Lyutikov Sz Blandford (2002), and the field in the early FS could be due to a small mixing 
of the ejecta with the shocked ISM. 

The transverse magnetic field in the ejecta decays as R~ l at early times, when the radial width 
of the material ejected in the explosion (SR) is constant, and the total energy in magnetic field 
is conserved. Therefore, if an explosion puts out equal amounts of energy in magnetic field and 
relativistic ejecta, the equipartition will continue to hold until 5R starts to increase. For a burst of 
duration T and Lorentz factor T, 5R ~ max{ cT, R/T 2 }, and so 5R increases only when i? cTT 2 . 
At large R, the magnetic field decays as R~ 2 and the energy in magnetic field decreases as 1/R. 
For GRB 021211, T ~ 4s & T ~ 500 (see §6), we expect a substantial fraction of the explosion 
energy in magnetic field at the deceleration radius of ~ 10 17 cm, if the burst was initially poynting 
flux dominated. 

We note that synchrotron radiation in the reverse shock could also explain the observed 7-ray 
fluence and the spectral peak frequency for 021211 (but not for the parameter space shown in fig. 2, 
which has too little flux at 50 kev). However, the low-energy spectral index in this case is ~ —0.5, 
whereas the observed index is 0.4 — f u oc v 0A — thereby killing this interesting possibility. 

The early and late time optical data is consistent with the density of the medium in the 
vicinity of the bust to be uniform, and a r~ 2 density profile is disallowed; the data can be fitted 
by a medium with 1/r 2 density profile provided that the density is smaller than normally expected 
for Wolf-Rayet stars at a distance of 10 cm by a factor ^ 10 4 . 

The upper limit on radio flux at 8.5 GHz at 10 days (Fox et al., 2003) poses a problem for 
the uniform density medium. The solution requires a combination of jet break at about 10 days 
and a decrease in es by a factor of ~ 2 between 11 min and 10 days, that we find not particularly 
appealing. 

We can set an upper limit of about 10 on the number of electron & positron per proton in the 
ejecta. The presence of pairs softens the reverse shock synchrotron radiation. When the charged 
lepton number exceeds about 10, the resulting reverse shock emission peaks at too low an energy, 
and produces flux at 8.5 GHz at 0.1 day that exceeds the observational upper limit by more than 
an order of magnitude. 

We thank Erin McMahon, Ramesh Narayan and Brad Schaefer for useful discussions. 
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Fig. 1. — The solid line is the ratio of the thermal energy per proton in the reverse shock and the forward 
shock, 7p,r/7p,/) as a function of n e j/no (the ratio of the comoving frame density of the unshocked ejecta and 
of the circumburst medium). The ejecta initial Lorentz factor, To, is used to normalize both these ratios such 
that the curves shown are independent of it. To a good approximation, To-y Ptr /j p j ~ 0.25(n e j /ncTo) -0 ' 7 • 
The relative velocity of the reverse shock front relative to the unshocked ejecta, as measured in the lab frame, 
V rs M b , is shown by the dotted line; r%V r3t iab ~ lA(n ej /Tln )~ - 5 . 
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Fig. 2. — Grey area: parameter space (and derived quantities) for a homogeneous external medium allowed 
by the observed R-band flux from reverse shock (RS) at 90 seconds and from forward shock (FS) at 11 
minutes, but not including the 7-ray flux during the burst. Black dots show the parameter space allowed for 
the observed gamma-ray flux to arise in forward shock as synchrotron emission; inverse Compton emission in 
FS or RS cannot account for the observed 7-ray spectrum. Top left panel shows the allowed density for the 
ISM (no) and e' e (which gives the minimum LF of shock heated electrons); gray area is allowed by the early 
reverse shock, and late forward shock, optical observations, and black dots show the region of the parameter 
space permitted by GRB observation. The top right panel shows the allowed value for the magnetic field 
parameter eg in the forward shock at 11 minutes. The lower left panel shows A v , a parameter that specifies 
the peak of the synchrotron spectrum at 11 minutes (the peak frequency is 4.95xl0 14 /A„ Hz). The lower 
right panel shows the ratio of es in RS and FS at 11 minutes (grey region), and the ratio of €3 in FS at 
deceleration time ( 3s) and at 11 minutes. Note that €b in RS is larger by a factor of about 10 3 compared 
with the value in FS at 11 minutes, and e_e in FS at deceleration is larger than at 11 minutes by a factor of 
about 10 2 . £52 = 30, td = 3s, p = 2.5 & z = 1.0 for all calculations; p — 2.2 gave very similar results. 



